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We study the spontaneous excitation of a radially polarized static multilevel atom outside a 
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6J0- 



spherically symmetric black hole in multi-polar interaction with quantum electromagnetic fluc- 
tuations in the Boulware, Unruh and Hartle-Hawking vacuum states. We hnd that spontaneous 
excitation does not occur in the Boulware vacuum, and, in contrast to the scalar held case, sponta- 
neous emission rate is not well-behaved at the event horizon as result of the blow-up of the proper 
acceleration of the static atom. However, spontaneous excitation can take place both in the Unruh 
and the Hartle-Hawking vacua as if there were thermal radiation from the black hole. Distinctive 

in ■ 

features in contrast to the scalar held case are the existence of a term proportional to the proper 

o" 

acceleration squared in the rate of change of the mean atomic energy in the Unruh and the Hartle- 
Hawking vacuums and the structural similarity in the spontaneous excitation rate between the 



cd 



^ static atoms outside a black hole and uniformly accelerated ones in a flat space with a reflecting 

H ■ 

boundary, which is particularly dramatic at the event horizon where a complete equivalence exists. 
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I. INTRODUCTION 

It is well-known that spontaneous emission (and excitation), as one of the most important 
features of atoms, can be attributed to vacuum fluctuations [l|, |2|], or radiation reaction [3j, 
or a combination of them [4j]. The ambiguity in theoretical interpretation, which roots 
in the freedom in the choice of ordering of commuting operators of the atom and field in 
a Heisenberg picture approach to the problem, was resolved by Dalibard, Dupont- Roc 
and Cohen- Tannoudji (DDC) (a, |6[ who proposed a formalism that demands a symmetric 
operator ordering so that the contributions of vacuum fluctuations and radiation reaction 
can be distinctively separated. The DDC proposal successfully resolves the problem of 
the stability of an inertial atom in its ground state in vacuum as a result of the delicate 
balance between the contributions of vacuum fluctuations and radiation reaction to the rate 
of change of the mean atomic energy. Recent investigations using the DDC formalism on the 
excitation of uniformly accelerated two-level atoms in interaction with fluctuating quantized 
massless scalar fields in vacuum in a flat spacetime with ml and without boundaries [8] and 
that of static atoms outside a Schwarzschild black hole [9] show that the delicate balance 
no longer exists in the cases under consideration, thus making the transition of the atom 
from ground state to excited states possible, i.e., excitation of atoms spontaneously occurs. 
The spontaneous excitation of uniformly accelerated atoms in the flat spacetime can be 



regarded as providing a physically appealing interpretation of the Unruh effect [lOj, since 
the spontaneous excitation of accelerated atoms gives a physically transparent illustration 
for why an accelerated detector clicks, while the spontaneous excitation of the static atoms 
outside a Schwarzschild black hole can be considered as providing another approach to 
the derivation of the Hawking radiation and it shows pleasing consistence of two different 
physical phenomena, the Hawking radiation and the spontaneous excitation of atoms, which 
are quite prominent in their own right. 

However, a two-level atom interacting with a scalar field is more or less a toy model, 
and a more realistic system would be a multilevel atom, a hydrogen atom, for instance, in 
interaction with a quantized electromagnetic field. Let us note that such a system in the 
multi-polar coupling scheme was recently examined in terms of the spontaneous excitation 



of an accelerated atom in both a free space U| and cavities 12|. It has been found that 



both the effects of vacuum fluctuations and radiation reaction on the atom are altered by 
the acceleration. This differs from the scalar field case where the contribution of radiation 
reaction is not changed by the acceleration. A dramatic feature is that the contribution of 
electromagnetic vacuum fluctuations to the spontaneous excitation rate contains an extra 
term proportional to a 2 , the proper acceleration squared, in addition to the usual Planckian 
thermal term of the Unruh temperature T = a/2ir. This is in contrast to the scalar field 
case where the effect of acceleration is purely thermal. As a further step along the line, in 
this paper, we would like to study, using the DDC formalism, the spontaneous excitation of 
a static multi-level atom in multi-polar interaction with quantized electromagnetic fields in 
vacuum outside a four dimensional Schwarzschild black hole. Our discussion will be based 



upon the Gupta-Bleuler quantization of free electromagnetic fields in a static spherically 
symmetric spacetime of arbitrary dimension in a modified Feynman gauge given by Crispino 



et al 13] . At this point, let us note the DDC formalism has also been app lied to calculate the 
radiative energy shifts of accelerated atoms both in a free space 14J4l6l| and in cavities J17F 
19 



The paper is organized as follows, we introduce, in Sec. II, the general formalism de- 



veloped in Ref. [8] and generalized in Ref. [11] to the case of a multilevel atom interacting 



with a quantized electromagnetic field in the multipolar coupling scheme. In Sec. Ill, we 
first review the Gupta-Bleuler quantization of free electromagnetic fields in the background 



Schwarzschild black hole [13], and we then define, in analogy to the scalar field case, the 
Boulware, Unruh and Hart le- Hawking vacuum states, calculate the two-point functions for 
the electromagnetic fields in these vacuum states and analyze their properties in asymptotic 
regions. The calculation of the spontaneous excitation of the multi-level atom interacting 
with a quantized electromagnetic field in the multipolar coupling scheme in all vacuum states 
will be performed in Sec. IV and a summary will be given in Sect. V 



II. GENERAL FORMALISM 

The DDC formalism is carried out in the Heisenberg picture that provides a very conve- 
nient theoretical framework as it leads, for the relevant dynamical variables, to equations of 
motion very similar to the corresponding classical ones. 

Consider a multilevel atom in interaction with vacuum electromagnetic fluctuations out- 
side a four-dimensional spherically symmetric black hole. The line element of the spacetime 

is 

d 2 s = (1 - 2M/r) dH - (1 - 2M/ry 1 d 2 r - r 2 (d6 2 + sin 2 6d(j) 2 ) (1) 

with M being the mass of the black hole. Let x(r) represent the stationary trajectory of 
the atom and r denote its proper time. The stationary trajectory condition guarantees the 
existence of stationary states. By assuming the multi-polar coupling, the total Hamiltonian 
that governs the evolution of the atom-field system with regard to the proper time r can be 
written as 

H(r) = H A (r) + H F (r) + H t {t) . (2) 

Here H^{t) is the Hamiltonian that determines the evolution of the atom, and it is given 
by 

h a{t) = y^a; n cr nn (r) , (3) 

n 

where \n) represents a complete set of stationary states of the atom with energies u n , and 
a nn — l n )( n l- Hp{r) is the Hamiltonian that decides the evolution of the free quantum 



electromagnetic field with respect to the proper time r, 

H F( T )=H UJ k a+ k a kJ^i ( 4 ) 

where k stands for the wave vector and polarization of the field modes, at and a^ are the 
annihilation and creation operators with momentum k. Hj(t) describes the coupling between 
the multilevel atom and the electromagnetic field. In the multipolar coupling scheme [201 ]. 

F 7 (r) = -e r(r) • E(x(r)) = -e J^ r mn • E(x(r)) a mn (r) , (5) 

mn 

where e is the electron electric charge, er the electric dipole moment of the atom, x{r) the 
atomic spacetime coordinates and E(a?(r)) the electric field operator of vacuum electromag- 
netic fields. 

Generally, atomic observables evolve with time as a result of interaction between the 
atom and the field. The rate of change of an arbitrary observable, 0(t), of the atom is 
governed by the Heisenberg equation 

^ = i[H A (r), 0(t)} - ie[r(r) • E(s(t)), 0(r)] . (6) 

Here we are interested in the second part on the right hand side of the above equation that 
is due to the interaction between the atom and the field 

^P) =-.e[r(r)-E(x(r)),0(r)]. (7) 

/ coupling 

The field operator E(a;(r)) can be divided into two parts as E = E-^ + E s . Here and after, 
the operators with superscript / represent the free parts that exist even when there is no 
coupling between the atom and the field and those with superscript s represent the source 
parts that are induced by the interaction between them. However, a tricky issue arises 
when we try to perform this decomposition in the above equation, since E-^ and E s do not 
separately commutate with the atomic observable. As a result, we can write 

^) = -ie{X&(x(T)) ■ ]r(r),0(r)] + (1 - A)[r(r), 0(r)] ■ W(x(r))) 

/ coupling 

-ze(\E s (x(r)).[r(r),0(r)] + (l-X)[r(r),0(r)}.E s (x(r))). (8) 

So, there exists an ambiguity of operator ordering. Dalibard, Dupont-Roc and Cohen- 

Tannoudji (DDC) [a, |6j proposed to use a symmetric ordering (A = |) so that the two terms 

on the right hand side can be separately Hermitian and possess an independent physical 

meaning that can be unambiguously identified as the contributions of vacuum fluctuations 

(vf) and radiation reaction(rr) respectively. Setting 0{t) to be the atomic energy H a {t), 

we have 

dH A {r) ( dH A {r) \ ( dH A (r) \ 

dr { dr ) vf + { dr J rr W 



with 



<IH {{T) ^ ^(Bf(x(T))-[r(r),^n<7nn(r)} + [v(T),J2^nn(r)}-Bf(x(r))) , (10) 



dr J . 2 

' V J n n 

d ^) = -f (E S (^W) • [r(r),^c n a nn (r)] + [r(r), ^ W „ ffnB (r)] • W(x(r))) . (11) 

Accordingly, we can also divide the dynamical variables of the atom, a mn or rj(r), into free 
and source parts, 

ff m»(T) = ^L( T ) + <m( r ) , ( 12 ) 

r l (r)=rf(r)+rf(r). (13) 

Then solving perturbatively the following Heisenberg equations of motion they satisfy 

-j-<Jmn(r) = i(u m - w n )a mn (r) - ieE(x(r)) ■ [r(r), cr mn (r)] , (14) 

— Ti(r) = i^2(u m - w„)(rj) mn cr mn (r) - ie[r(r),ri(r)] • E(x(r)) , (15) 

mn 

to order e, we find 

u mn\ l ) u mn\'0) c i ttn\ 

(16) 

<n(r) = ^ej; o dr'Ef(x(r'))-[rf(r'),aUr)}, 

(rf (r)) mn = (rf (r )) m „e 1 ^— )(— ") , 
r?(r) = -.e r;rfr'[r|(r'),rf(r)]Ej(x(r')) . 

Here the repeated subscript denotes the summation over all spatial components and 
(rf(r)) mn = (m|r^(r)|n) . Similarly, using the Heisenberg equation for the annihilation 
operator of the field, 

^ a Mr)) = -iuj n a n (t(r))-ier(r) ■ [E(x(r)),o £ (t(r))]^ . (18) 



one can show to the same order of perturbation that 






4(*( r )) = 4W r °)) e ~ jus 



al(t(r)) = -iej T n) dr' r-f(r') • [W(x(t% a^t(r))} 



(19) 



and therefore the source filed E s that is generated by interaction between the atom and the 
free field can be expressed as 

E?(s(r)) = -ie /VrJ(r')[Ej(:r(T')),Ef(:r(T))] . (20) 

J TO 



Now we suppose that the field is a vacuum state and the atom is in state |6). For 
simplicity, we also assume that the atom is polarized along the radial direction defined by 
the position of the atom relative to the black hole space-time rotational killing fields. This 
assumption significantly simplifies, while keeping the physical essence of the problem, the 
computations we are going to perform since we do not need to calculate the contributions 
associated with the polarizations in the 9— and 0— directions. Taking the expectation values 
of Eqs. ( TTU1) and (TTT1) over the state of the system |0, b), we can, using Eqs. ( JT6l) . ( TT71) and 
(T20J) . show that the contributions of vacuum fluctuations and radiation reaction to the rate 
of change of the mean atomic energy are given, to order e 2 , by 

l 6 ^-) = 2ie 2 £dr' C F (x(T),x(r'))± X t(T,T') , (21) 

(^^) =2.e 2 £rfr' X ^(x(r),x(r'))^(r,r'). (22) 

In the above equations, C F and \ F are the two statistical functions of the electromagnetic 
field, i.e., the symmetric correlation function and the linear susceptibility function. The 
radial components of these function, which are relevant in our future calculations, are defined 

as 

C F (x(r),x(r')) = i(0|{£/(*(r)),£/(*(r'))}|0> , (23) 

X F (x(r),x(r')) = ±(0\[E}(x(T)),EfaW)]\0) > (24) 

where { , } denotes the anti-commutator and |0) represents the vacuum state of the field 
which will be defined in the next Sect. Two statistical functions of the field are dependent 
on the trajectory of the atom. Analogously, C^(r, r') and xt( T ^ T ') are two atomic statistical 
functions which are defined as 

C b A (r,r') = ^|{r'(r),rV)m, (25) 

Xt(r,r') = l(b\[v f (r),rf(r')]\b). (26) 

They don't depend on the atomic trajectory and are determined only by the internal struc- 
ture of the atom itself. Their explicit forms are given, with respect to r, by 

C*(t, r') = -J2\ ( & l r (°) \ d ) I 2 [e iu)bd{T - T,) + e -^( T - T ')] , (27) 

d 

Xb(r,r') = ^K 6 l r (°)l rf >| 2 [^ M(T - r,) - e-* 1 *^)] , (28) 

d 

here w« = u\, — Ud and the sum extends over the complete set of atomic states. 

Now it is clear that the calculation of the rate of change of the mean atomic energy 
requires detailed knowledge on the quantization of electromagnetic fields in the exterior 
region of the black hole and the specification of vacuum states. This is the main topic for 
the next section. 



III. QUANTIZATION OF ELECTROMAGNETIC FIELDS IN THE EXTERIOR 
REGION OF A SCHWARZSCHILD BLACK HOLE 



The quantization of the electromagnetic field in an static spherically symmetric 
Schwarzschild-like spacetime has been carried out by Crispino et al [13| using the Gupta- 
Bleuler condition in a modified Feynmann gauge. Here we first give a brief review of their 
basic results, and we then define the vacuum states and calculate the statistical functions of 
the field. The Lagrangian density of the electromagnetic field in a modified Feynman gauge 
is 



£, 



-g 



4^ r 



2 G 



(29) 



with G = VM M + K^A^ and K^ being a vector independent of A^. From the Lagrangian 
density, we can write down the field equations. If we choose K^ to be K^ = (0, 2M/r 2 , 0, 0), 
the equation for A t decouples from other ones 13] . A complete set of solutions of the field 



equations can then be denoted by A^ n,u) . Here the label "n" distinguishes between modes 
incoming from the past null infinity J~ (denoted with n =+- ) and those going out from the 
past horizon H~ (denoted with n =—>■)■ The modes with A = are 



a (On; culm) 



R? n \u\r)Y, 



lm ^ 



-iuit 



0,0,0) 



where Y$ m is the spherical harmonics and R t n (u\r) satisfies the radial equation 
■- l l - 2M/r) d ( nd\ 1(1 + 1) 



uj- 



2M/r] 



+ 



dr \ dr 



R(° n )(a;|r) = 0. 



(30) 



(31) 



These modes are nonphysical as they do not satisfy the gauge condition G = that are 
satisfied by all other modes with A 7^ 0. Modes with A = 3 are pure-gauge modes, and they 
are given by 

A(3n;u)lm) vy \{nuilm) ( r { r >\ 

with 

k {nujlm) = - Rf n \u\r) Y lm e~ iujt . (33) 

to 

Those with A = 1,2 correspond to two classes of physical modes. For the first class of 

physical modes (A = 1), A t — 0, and 

A {iwMm) = j R( 1 ")( w | r ) y lm e- twt 
with / > 1, where the radial function R t n (co\r) satisfies 

>lln) (u;|r) = 0. (35) 



(34) 



r 2 dr 



d 



„2 D (ln), 



(l-2M/r)j-(r 1 RY n \uj\r)) 



+ 



u> 



1(1 + 1) 



2M/r) 



R 



The angular components can be expressed as 



.(In; o,Im) _ l~2M/r d f r 2 R (ln) (l) \ 8Y 

A { - ^ dr \rR t [u\r)jd,Y lm e 



-iuit 



(36) 



For the second class of physical modes (A = 2), A t = A r = 0, 

Af n ' ulm) = Ri 2n \co\r) if m) e' iwt (37) 

with R\ n '(uj\r) obeying the following radial equation 



co 2 d ( d\ 1(1 + 1) 



+ — (1 - 2M/r 



1 — 2M/r) dr \ dr J r 



.2 



>(2n) 



i?P(w|r) = (38) 



in which Y { are the divergence-free vector spherical harmonics on the unit 2-sphere sat- 
isfying 

V k (V k Y} lm) -V t Y^ m) ) = -l(l + l)Y l {lm) , (39) 



dQ fj ij r/ /m) if m ' } = s lv s mm/ . (40) 

Here the overline denotes the complex conjugation, % angular variables on the unit 2-sphere 
S 2 with metric fjij and inverse metric ffi with signature (+, +), Vj the associated covariant 
derivative on S 2 and V* = fj^Vj. The above four classes of modes form a complete set 
of basis for the quantum electromagnetic field. The normalization of them are determined 
from the canonical commutation relations of the fields by requiring suitable commutation 
relations for the annihilation and creation operators. 

To quantize the electromagnetic field, let us define the general inner product as 

(A®, A^) = f d^ W[A^\ A«)\ (41) 

in which 



w*[A«\ a^')} = i [ a { P nK')"" - n(o^ aP ] , (42) 

w ^7=iwk] = - [F, "' + ' rG] (43) 

and dYin = da n M , Aj2 = A^ n,0J . The equal-time commutation relations for the fields and 
their momentum operators are 

[ i M (t,x), i„(t,x') ] = [ n^(t,x),n^(t,xO ] = , (44) 

[ i M (t,x),n^(t,x') ] = -j£= o 3 (x- x') , (45) 

where x and x' represent all spatial coordinates. Expand the field operator in terms of the 
complete set of basic modes as 

A,(t, x) = W°° -^ l4 Xn ^ lm Ht, x) aS + A^ lm \t, x) a + Jt n) \ , (46) 



Xnlm 



By using the inner products of the field and the commutation relations between the field 
and momentum operators, the commutation relations between the annihilation and creation 
operators are found to be 

[«Sm> «L 3 '/w] = -[ a Sm> &m'\ = <WMmm'<*(w - £*/) , (47) 

[«2m> a S"w] = [«Sm> «2"'m'] = <W<^<W'<*(u; - w') (48) 



and all other commutators vanish. The Gupta-Bleuler condition [21] requires that 

G+|p%s) = 0, (49) 

where G + is the positive-frequency part of G — VM M + K^A U and \phys) represents an 
arbitrary physical state. This is equivalent to the following statements: the states obtained 
by applying the creation operator aj^l are all nonphysical; the physical states of the form 
<K)frn \phys) have zero norm and are orthogonal to any physical state. 

Now we can define vacuum states and calculate the correlation functions of the field. We 
start with the Boulware vacuum state, |0)b, which is defined by 

a Lm|0)-B = j f or a tt (A, n, u;,Z, m) with (a> > 0) . (50) 

In the computation of the correlation function B(0\E r (t)E r (t')\0) b with E r = A t - r — A r]t , 
the contributions of nonphysical modes and pure-gauge modes are found to be canceled out, 
and only the contribution of the first class of physical modes is left 

\E r (x)E r (x')\0) B = —J2 / dtuue^-^ Y lm {6 lV ) Y* m (9',v') 

[^,(w|r)^?(w|r / ) + S,(w|r)^ r (w|r')] . (51) 



lm 

XI 



Here and after, for the sake of brevity, we omit the label A = 1 as others will not appear in 
the correlation functions. For the Unruh vacuum state, \0)u, similar to that in the scalar 



field case [1(J, |22|, we define the positive frequency for modes steming from H~, the past 
horizon of the black hole, with respect to the Killing vector £ = d\j and those originating 
at infinity with respect to the Killing vector rj = d t . For the Hartle-Hawking vacuum state, 
\0)h, we define the incoming modes to be positive frequency with respect to v, the canonical 
affine parameter on the future horizon, and the outgoing modes to be positive frequency with 
respect to u, the canonical affine parameter on the past horizon. Repeat the same steps as 
that in the Boulware vacuum case, the correlation function satisfying the corresponding 
boundary conditions in both the Unruh and Hartle-Hawking vacuums are also found to be 
only associated with the first class of physical modes. They are given respectively by 

\E r (x)E r (x')\0)u = —J2 / duuje-^ t -^Y lm (e,<p)Y l U6\<p') 

'^,(w|r)^f(w|r / ) 



l in 



X 



-2ttui/k 



+ 9(cu) < R l (co\r) < Rt(^\r') 



(52) 



and 

H (0\E r (x)E r (x')\0) H 



-E 

471- A^ 
Im 



du U> 



e y /m (#,</)) Y lm \p ,ip )- 



1 



-2-klu/k 



+ e ^ m (0, ¥?J r, m ^ , v J — e 2^/ K _ i — 



3' ,_^-"*' 



<=r 



(53) 



where k = 1/4M is the surface gravity of the black hole. Details about the features of the 
correlation functions calls for the specific properties of the radial functions, so now we turn 
our attention to the analysis of the radial functions. To do so, let us further write the radial 
function as 



>(«) 



*T>|r) 



v^+iJv&V) 



CO 



then Eq. (135]) becomes 



^ 2 
drl 



2M\ 1(1 + 1) 



» 



<P%\r) = 0, 



(54) 



(55) 



ir 1 / iviii 

where r* = r + 2Mln(r/2M — 1) is the Regge- Wheeler tortoise coordinate. It is difficult to 
exactly solve this equation in terms of elementary functions, but fortunately, two classes of 
physical solutions in the two asymptotic regions single out 



~$ui{r) 



e iujr " + J?i(u) e 



■^r 



<Pui(r) 



Ti{uj) e* 
Ti(uj) e~ 

e -iur, + 



%{uj) 



r^2M 
r — >• oo , 

r^2M 
r — i- oo . 



(56) 



(57) 



Here 1Z and T are, respectively, the reflection and transmission coefficients and the following 
relationships exist among them and their complex conjugates, 



fi(") 



^M 



(58) 



Ti(u), 

$i(u)\ = \fci(u)\ , 

1 - &(u)\ 2 = 1 - \%(u)\ 2 = \T t (u)\ 2 , 

k rfi*(u)Ti{«>) = -Tf{uj)%{uj) . 

Although the exact forms for these coefficients demand an exact solution of Eq. (153]) which 
is a formidable task, further studies (see Appendix) show that the summation concerning 
the radial functions in the two asymptotic regions, i.e., r — > 2M and r —¥ oo, behaves as 



i- 



E, *('+l)(2l+l) |7I(a 



and 



^2(21 + 1) \R t (uj\r 



^,{2l + l)\%(u)\r) 



(2M) 4 u> 2 



3 9(M 



8c/ 



+ 



1 

&M 2 g'L ' 



3g 00 ™j aoo 
E,t(i+l)(2t+l)|7IM| 2 



r ^2M 



r — > oo 



r ^2M 



r — ► oo . 



(59) 



(60) 
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IV. SPONTANEOUS EXCITATION OF AN ATOM INTERACTING WITH VAC- 
UUM ELECTROMAGNETIC FIELD IN SCHWARZSCHILD SPACETIME. 

Suppose that a static multilevel atom is held at the radial distance r and it interacts with 
fluctuating quantum electromagnetic fields in vacuum. Using the DDC formalism introduced 
in section II and the results in the preceding section, we now calculate the rate of change of 
the mean atomic energy. 

a. Boulware vacuum. In the Boulware vacuum state, the two statistical functions of 
the electromagnetic field, i.e., the symmetric correlation function and the linear susceptibility 
function, can be easily found from Eq. f[5"Tj) to be 



C (x(t),x(t')) = -^-^ / du u (e v^oo+eV^oo) 

'o 

00 



327T 2 



^(2l + l)[\%(u\r)\ 2 + \%(u\r)\ 2 ], (61) 



x 

1=1 



and 



iuj/S.T zcjAt 



X (x{t),x{t')) = -— / duu{e vm-evm' 

OZ7T Jo 



oo 
X 
1=1 



^(2l + l)[\%(u\r)\ 2 + \%(u\r)\ a ]. (62) 



In obtaining the above results, we have used the relation At = y/g~oo At with g 00 = (I — 
2M/r) and the following property of the spherical harmonics 

El^(^)! 2 = ^- ( 63 ) 

m=— I 

Inserting Eqs. flBTT) and (T2"5j) into Eq. ( l2~Tj) . we obtain the contribution of vacuum electro- 
magnetic fluctuations to the rate of change of the mean atomic energy, 

' vf \- Ub>UJd 

- E \(b\r(0)\d)\ 2 u 2 bd P(-u; bd ,r) 

where we have defined 

P{u, r) = ~P(co, r) + P~(cu, r) (65) 

with 



(64) 



f(u, r) = ^(2* + 1) \t(uy/g^ \r)\ 2 , (66) 

i 



11 



Hereafter, the summation over / is implied to range from 1 to oo. Similarly, by inserting 
Eqs. ( E2J) and (12"T|) into Eq. ( I2"2"j) . the contribution of radiation reaction to the rate of change 
of the mean atomic energy is calculated out to be 

•dH A (r)\ «"*•[£ |<Hr(0)|<OI'<4n**r) 



dr /„„ 167r 

+ £ K&|r(0)|d>| 2 u&P(-u; M ,r) 



(68) 



Adding up Eqs. (1M|) and (I6"8"j) yields the total rate of change of the mean atomic energy 
'dH A (r)\ e 2 g 00 



£ l(&|r(0)M)| 2 a; 6 2 d P(a; M ,r). (69) 



dr / .. 8ir 

101 U b >UJ d 

Here only the term with ui^ > u^ survives after the addition and it is negative. This means 
that for an atom originally in the ground state, the contributions of vacuum fluctuations 
and radiation reaction to the total rate of change of the mean atomic energy cancel, and, 
as a result, the ground state is stable. However, an atom originally in an excited state can 
transition to lower level states, since the total rate of change of the mean atomic energy 
is negative. Although, qualitatively, these features are same as those for the atom in the 
Minkowski vacuum state in a flat spacetime, the total rate of change also displays quanti- 
tative difference which is embodied in the factor P(ojbd,r). A comparison of Eq. (|69|) with 
Eq. (23) in Ref. [7|, which gives the rate of change of the mean atomic energy for an inertial 
atom in a flat space with a reflecting boundary, shows that the two rates are quite similar, 
and the appearance of P(ujm, t) in Eq. (169]) can be understood as a result of backscattering 
of the vacuum electromagnetic field modes off the space-time curvature in much the same 
way as the reflection of the field modes at the reflecting boundary in a flat space-time. To 
further understand the effect of backscattering caused by the spacetime curvature, let us 
now analyze what occurs in two asymptotic regions, i.e., at the spatial infinity and at the 
event horizon, which are regions of particular physical interest to us. 

At spatial infinity, i.e., r — > oo, a further simplification of the total rate of change by use 
of Eqs. ([59]) and @0J yields 

' tot w b >w d 

in which _ 

,, , 3E^(^ + 1)(2/ + 1)|77(^ V ^)| 2 

/("«> r ) = ^Ua ( 71 ) 

is a grey-body factor that characterizes the backscattering of the electromagnetic field modes 
off the spacetime curvature. Notice that at infinity, /(ww,r) ~ 0, so the rate of change re- 
duces to that of an inertial atom in the Minkowski vacuum in flat spacetimes with no bound- 
aries, suggesting that the Boulware vacuum at large radii is equivalent to the Minkowski 
vacuum. 
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When the atom is fixed near the event horizon, i.e., when r — > 2M, further simplification 



gives 



with 



dH A (r) 

(It 



tol 



e 



E [<6K0)|d>| 



2 4 



aj b >aj d 



M 



bd 



M 



1 + 



W 



/;-/ 



+ /(ww, r) 



(72) 



(73) 



/fl-oo ryi -2M/r ' 

where a is the proper acceleration of the static atom. Notice the appearance of an extra 
term proportional to a 2 here in contrast to the scalar field case [9[. Let us note however 
that the presence of the terms proportional to the proper acceleration squared also appear 
in cases when the Unruh-DeWitt monopole detector is replaced by a dipole detector which 
couples to the derivatives of a scalar field [23|, |24J. The proper acceleration a diverges as the 
event horizon is approached, so does the rate of change of the mean atomic energy. However, 
at infinity, where the spacetime is asymptotically flat, a ~ 0, so its contribution to the rate 
is negligible. This result is in sharp contrast to that of the scalar field case, where the rate 
of change is always finite [9(. Except for the grey-body factor, /(w M ,r), here the result 
agrees with that of a coaccelerated atom with a proper acceleration a in interaction with 
fluctuating electromagnetic fields in the Rindler vacuum 25] (the case with the temperature 
of the thermal bath set to zero). The above discussions reveal that the Boulware vacuum 
is the vacuum state of static observers outside a black hole, and it resembles the Rindler 
vacuum in the flat spacetime since the static atom is accelerating with the proper acceleration 
a with respect to locally free-falling observers. 

b. Unruh vacuum. For a multilevel atom in interaction with quantum electromagnetic 
fluctuations in the Unruh vacuum, two statistical functions of the field are easily obtained 
from Eq. ( 152]) as follows 



C F (x(r),x(r')) 



1 



327T 2 



x 



juj/S.t jujAt 

dw u I e v'soo + e v^oo 



E,(2/ + l)|^Hr-)| 2 

1 p — 2nuj/K 



6(u)^2(2l+l) \%{w\r 



(74) 



and 



X F (x(r),x(r')) 



327T 2 



zjajAt ico At 

du uj I e v'soo — e v^oo 



E,(2/ + l)|^Hr)| 2 

I p—2ttlu/k 



9(cu)J2( 2l + 1 ) I^M r 



(75) 



Now the contributions of vacuum fluctuations and radiation reaction to the rate of change of 
the mean atomic energy can be calculated by inserting the statistical functions into Eqs. ( 12~T|) 



13 



2 , .2 



and (|22|) . For the contribution of vacuum electromagnetic fluctuations, we have 

dH A (T)\ 

I vf 

1 W, P(-U bd ,l) ■(„ 



dr 



e 2 goo 
' 16tt 



E l(6l r (0)|d>| 



^6d X 



1 + - 



2KL} bd /K r 



~^(ubd,r) + ;' J"*' ' ' -- ^(w 6d ,r) 



e 27r ^6d/ K r _ x 



|2 ^x 



- E K*i r (°)i<oi : 

1 N-^/ s , P(fcW 



^6<W d 
1 + 



e 27r|o; 6d |/K r . 



2w\ui bd \/K r 



+ P(-u bd ,r) 



where we have defined 



y/l - 2M/r ' 
For the contribution of the radiation reaction, 

dH A {r)\ e 2 g 00 f 

/rr ^ 167r I 

rT K LU b >LU d 

1 + - 



(76) 
(77) 



dr 



E i^i r (°)i 



1 2 , .2 



W M X 



D 2lTbJ bd /K r _ ]_ 



?<■ 



w M , r ) 



P(-w 6d ,r) <- 

- + P {uj m , r) 



D 2lTL0 bd /K r 



12 <4d x 



E l<6|r(0)|d>| : 

^— - j P (-ww, r) - e27r| , M|/Kr _ x + P (-ww, r) j | • 



Ul b <U! d 
1 + 



3 27t|o; m |/k 



(7* 



Compared with the case of the Boulware vacuum, both the contributions of vacuum electro- 
magnetic fluctuations and radiation reaction are altered by the appearance of a Planckian 
factor. Adding them up, we obtain the total rate of change of the mean atomic energy, 



dHAr) 

dr 



2 9oof 



tot 



8/r - E \mo)\d)\ 



2 2 
U bd 



• E I WW I 

u b <u d 



2 2 

CO 



g27ro; M /K r — X 

P{uJ bd ,r) 



P(u bd ,r) + P(u bd ,r) 



bd 2ir\u bd \/Kr 



(79) 



Now the delicate balance between vacuum fluctuations and radiation reaction that ensures 
the stability of the atom in its ground state in the Boulware vacuum no longer exists. The 
u b < bj d term which gives a positive contribution to the total rate of change of the mean 
atomic energy makes the transition of the atom from the ground state to an excited state 
possible, i.e., excitation spontaneously occurs in the Unruh vacuum outside a black hole. The 
structure of the rate of change also suggests that there is thermal radiation from the black 
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hole (represented by the Planckian term) which is backscattered by spacetime curvature 
(represented by r). It is this thermal radiation that renders the spontaneous excitation 
possible. The temperature of the thermal radiation is given by 

T = £ = £ /i \m, = ^r 1/2 Tn , (80) 

Itx Lix y/i - 2M/r 

with Tn — k/27t being the usual Hawking temperature of the black hole. This is actually 
the Tolman relation which gives the temperature felt by a local observer. 

In order to gain more understanding, let us now examine the behavior of the rate of 
change in two asymptotic regions. First, when the atom is fixed near the event horizon, i.e., 
when r — > 2M, we can approximate, by use of Eqs. fl59l) and (160]) . the rate of change as 
follows 



cIHa(t)\ e 



dr /,. 37r . 

101 K U) b >U) d 



|2 , ,4 
W bd 



l + : Jt 1 H 1 + :^)+/Uv,,/.r.) 

bd, 



e Wbd l T — 1 / V w? ■■ 



1, *•,,> , \ bd / J 



A distinct feature in contrast to the case of the scalar fields [9| is the existence of an extra 
term proportional to a 2 , the proper acceleration squared. It is worth pointing out here that 
the appearance of such a term has also been found when one studies the spontaneous excita- 
tion of an accelerated multilevel atom coupled with electromagnetic vacuum fluctuations in 



a flat spacetime [111 ]. Noteworthily, one can show that, close to the event horizon, T rs a/2ir 
holds. This leads to a remarkable observation, that is, the u^ < u>d term, which makes the 
spontaneous excitation possible, can be viewed completely as a result of the Unruh effect 
due to the proper acceleration that must exist to hold the atom static (refer to Eq. (29) in 



Ref . [Uj ) . This demonstrates a complete equivalence between the effect of acceleration and 
that of a gravitational field at the event horizon. 

If the atom is placed far away from the black hole, i.e., r — > 00, the total rate of change 
of the mean atomic energy becomes 

'd H A(r)\ _ e 2 f ^ u ,,, n .,«,2.A 



dr I ,. 3ir ■ 

101 v iJi,>OJ c J 



J2 l<6|r(0)|d)|-a; M 



1+ ( 1 + ^w^— t )/(<""> ') 



-EIWMI'^ ^U . (82) 

U) b <U) d 

Notice that thermal terms are all multiplied by a grey-body factor, /(u^, r), which vanishes 
at spatial infinity. The above results are in accordance with the common belief that thermal 
flux emanates from the black hole horizon and is backscattered and partly depleted by the 
curved spacetime geometry on its way to infinity. So, as the atom is placed further and 
further away, the flux it feels becomes weaker and weaker. 
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c. Hartle- Hawking vacuum. Now let us look at the case in which the electromagnetic 
fields are in the Hartle-Hawking vacuum state. Two statistical functions of the field can be 
easily found from Eq. ( 153]) to be 



C f (x(t),x(t>)) 



32tT 



iojAt zljAt 



duo uj e v^oo + e v^oo 



X 



£,(2J + 1) \&Mr)\l + £,(21 + 1) \%(u\r)\ 2 



1 o—2ttuj/k 



■,2-ttlu/k 1 



(83) 



and 



X F (x(t),x(t')) 



32tt 2 



zuAt jujAt 



dcu uj [ e ^m _ e vw 



x 



£,(2/+l)|if,Mr)| 2 ^(2/ + l)|^(^k)| 2 



1 p—2-rru/K 



2ttlj/k 1 



(84) 



Inserting them into Eqs. (|2T|) and (122]) yields the contributions of vacuum fluctuations and 
radiation reaction to the rate of change of the mean atomic energy 



and 



dH A (r) 
dr 



vf 



e 2 9oo 
16tt 



E K 6 l r (°)l d >l 



2 2 
U bd 



UJ b >Ul d 



P(-u bd ,r) + P{u bd , 



g2iru) bd /n r — J 



+ 



1 + 



^2-wbj bd /n r — 1 / V 






(^(w M ,r) + £(- 



■wmi r 



E |<6|r(0)|d)| 2 u; 6 2 d 

1 

1 + 



w 6 <aj d 



?< 



t^, 



ww,r) + P(-o; M ,r) 



+ 



e 27r|^bd|/K r _ ]_ 



e 27r|w M |/«r _ ]_ 

~?(-Ubd,r) + P(ubd,r] 



/ dH A (r) 
\ dr 



e 2 9oo / 



167T 



E K6|r(0)|d>| 

V/.J, S./.1 , 



2 , , 2 
w bd 



-UJ b >Ul d 



P{-u bd ,r) + P{u bd ,r) 

p2-nui hd /n r — -y 



1 + 



^nw bd / 



j^—A(^{u bd ,r) + P~{-u M ,r)\ 



2 2 

tot 



+ e \(b\m\d)\ - w 
i 



oj(,<o; d 



^(^&d, 



r) + P(-u bd ,r) 



1 + 



j27r|aj bd |/K r _ ]_ 



e 27r|oj ftd |/K r _ ]_ 

P(-w M ,r) + P(w M ,: 



55) 



(86) 
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Then total rate of change readily follows 



/ dH A (r) 

\ dr i ,„, 



H E iw°mi 



8tt \ 



2 ^L x 



UJ b >UJ d 



1 + 



-i_ I j^( WM> r) + ^(-w W> r)J 

E |(6|r(0)|d)| 2 a; 6 2 d x 



w b <w d 



o2wui bd /n 



(87) 



As in the Unruh vacuum case, the positive uj\, < ujd term also appears, and this term leads to 
spontaneous excitation of static atoms in the Hartle-Hawking vacuum in the exterior region 
of the black hole. Besides, the Planckian factor in the total rate of change is a revelation 
of the thermal nature of the Hartle-Hawking vacuum. To learn more, let us further study 
what happens in the asymptotic regions. 

When the atom is fixed at infinity, i.e., r — y oo, the rate of change of the mean atomic 
energy is 



dHAj) 

dr 



tot 



3tt 



f E iw°)i rf )i 2 4 1+ 



LU b >bJ d 



1 



■ E K & l r ( MI 



2 a; 4 



g27rai M /K r . 
1 



l + /(ww,r) 



bd 2-K\w bd \JK r 



UJ b <UJ d 



1 + /(«m, r) 



(88) 



furthermore, at spatial infinity where f(ubd, r ) can be taken as zero, the total rate of change 
is what one would obtain when the atom is immersed in a thermal bath at the Hawking 
temperature T = Tjj in a flat spacetime, while at the event horizon, i.e., r — y 2M, 



/dHAj) 
\ dr 



i,, i 



f E K*|r(0)|d>r<4i(l + 



UJ b >U) d 

■E 



K&K0)|rf)| 2 ^ 



e 27r ^6d/ K r 
1 



w 



/;(/ 



1 + ^ +/(WM)0 



^Trl^dl/K,. 



W, 



/;</ 



1 + ^ +/(wm>0 



,(89) 



it is divergent as a — )■ oo. In addition to the contribution of the outgoing thermal radiation 
from the event horizon that exists in the Unruh vacuum (refer to Eq. flHTj) ). there is another 
contribution to the total rate of change, the thermal term multiplied by f(ubd, r ), which can 
be regarded as resulting from the incoming thermal radiation from infinity. Both incoming 
and outgoing thermal radiation are backscattered off the spacetime curvature on their way. 
This result is consistent with our usual understanding that the Hartle-Hawking vacuum is 
not empty at infinity but corresponds instead to a thermal distribution of quanta at the 
Hawking temperature, thus it describes a black hole in equilibrium with an infinite sea of 
black-body radiation. 
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A distinctive feature in contrast to the scalar field case is the existence of the term 
proportional to a 2 that is nontrivial at the event horizon and thus is physically important. 
Remarkably, Eq. f|89|) is in structural similarity to the total rate of change of the mean atomic 
energy of an uniformly accelerated atom interacting with electromagnetic field fluctuations 



in a flat space with a reflecting boundary (refer to Eq. (60)) 12]), reflecting again that the 
scattering of the electromagnetic field modes off the spacetime curvature plays similar role 
as the reflection of the field modes at boundaries in a flat spacetime. This similarity is 
particularly striking at the event horizon, where T ~ a/27r. 



V. SUMMARY 

In summary, using the Gupta-Bleuler quantization of free electromagnetic fields in a static 
spherically symmetric spacetime of arbitrary dimension in a modified Feynman gauge given 



by Crispino et al [13], we have defined, in analogy to the scalar field case, the Boulware, 
Unruh and Hartle-Hawking vacuum states outside a four-dimensional Schwarzschild black 
hole, calculated the two-point functions for the electromagnetic fields in these vacuum states 
and analyzed their properties in asymptotic regions. We then computed the contributions 
of vacuum fluctuations and radiation reaction to total rate of change of the mean energy 
for a radially polarized static multilevel atom in interaction with quantum electromagnetic 
fluctuations in all the three vacuum states. 

Our results show that the static atoms in the ground state in the Boulware vacuum are 
stable and this is in qualitative agreement with the case where the atom is assumed to be in 
interaction with quantized massless scalar fields [9]. However, the spontaneous emission rate 
of the excited atoms close to the horizon contains an extra term proportional to the squared 
proper acceleration of the atom in contrast to the scalar field case, and this rate is not 
well-behaved at the event horizon as a result of the blow-up of the proper acceleration of the 
static atom with respect to the free-falling local observers there (note that this acceleration 
vanishes however at the spatial infinity). This is in sharp contrast to that of the scalar field 
case, where the rate of change of the mean atomic energy is always finite [9J. 

For the static atoms in both the Unruh and the Hartle-Hawking vacuums, the delicate 
balance between the contributions of vacuum fluctuations and radiation reaction that ensures 
the stability of the static atoms in ground state in the Boulware vacuum no longer exists, 
so spontaneous excitation occurs in the exterior region of the black hole. The spontaneous 
excitation rate of the static atoms is in accordance with our usual understanding that the 
Unruh vacuum describes a black hole with thermal radiation emitting from its event horizon, 
whereas the Hartle-Hawking vacuum depicts a radiating black hole in equilibrium with an 
infinite sea of black-body radiation. 

Distinctive features in contrast to the scalar field case are the existence of the term 
proportional to the proper acceleration squared in the rate of change of the mean atomic 
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energy in the Unruh and the Hartle-Hawking vacuums and the structural similarity in the 
spontaneous excitation rate between the static atoms outside a black hole and uniformly 
accelerated atoms interacting with electromagnetic field fluctuations in a flat space with a 
reflecting boundary, which is particularly dramatic at the event horizon where a complete 
equivalence exists. 
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Appendix: on the asymptotic evaluation of mode sums 

Here, we derive the properties of the quantity 5^(2/ + 1) \R\ , (tu\r)\ 2 in two asymptotic 
regions. Let us start with the incoming modes Ri(u\r), by examining, at a fixed radial 
distance r, the correlation function of the field in the Boulware vacuum, Eq. ( 15TT) 



1 POO ^—_, i. 

B (0\E r (r)E r (r')\0) B = — - / du to V(2Z + 1) [ \R l 



[co\r)\ 2 + \Ri(u\r)\ 2 } e v^oo 



' (A.l) 



7T 2 (At) 



Here At is the interval between the coordinate time and the approximation is taken at 
spatial infinity, where the spacetime is asymptotically flat. So, the two point function in the 
Boulware vacuum at large radii and that in the Minkowski vacuum should agree. At spatial 
infinity (r — > oo), it can be deduced from Eq. (156]) that 

B^ + l)I^M0l^ E '' ( ' + 1)(2 ^ 1)|7iM|2 . (A.2) 



u'r H 



This is very small as r — > oo, and thus can be neglected in Eq. (1A.1J) . Using f u 3 e 



-tux b 



we obtain 

,2 



£(2Z + l)|^|r)| 2 ^!^r. (A.3) 
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The summation over the outgoing modes, Ri(co\r), in the region r ~ 2M can be obtained 
by solving the equation of the corresponding radial function as follows. Let £ 2 = r/2M — 1 
and q = 4Mu, the radial equation Eq. (153]) can be simplified to be 

e d ^+^ + [ q 2 -(2lO 2 tfi = 0, (A.4) 

we have approximated 1(1 + 1)£ 2 by (Z£) 2 since £ ~ 0. The general solution of this equation 
can be expressed in terms of the modified Bessel functions as 

^i\r-,2M ~ a t K iq (2l£) + k 7_ <g (2/0 • (A.5) 

To estimate the coefficients a\ and bi, let us look at the radial equation Eq. (1551) in the 
limiting case, I — >• oo. Now the effective potential in the equation is very large as compared 
to the other two terms for fixed r and u, therefore, one can deduce that ipi ~ for large I. 
As a result, bi is an exponentially small function of I when I is large, as J_j g (2/£) ~ f 
(/ 3> 1). The second part in Eq. flA.5[) therefore makes a bounded contribution to the 
summation of 5^(2/ + 1) | Iii(uj\r)\ 2 and it is negligible as the contribution of the first term 
is proportional to £ -4 . The coefficient a\ can then be determined by comparing the general 
solution, Eq. (IA.5J) . with the asymptotic solution 

^(r)~e"'+4)e-"', (A.6) 

and the result is 

2 t-iq e iq/2 

Thus the summation to the leading order is 

^(2/+i)|itMr)i 2 ~ ^ r4r(2g 4 )rHg) E^+ 1 )( 2/ + 1 )i^™r 

8a; 2 1 

W 3^ 2 + 6M 2 ^ 2 • ( • j 

Here we have appealed to the trick of approximating the infinite summation over I by an 
integral. 

The summation over the outgoing modes, Ri(tu\r), at infinity can be easily deduced from 
the asymptotic solution of the outgoing mode, Eq. (15BJ1 as 

B2i + Dl^k)l^ E '' ( ' + 1)(2 2 ': 1)|7;M|2 . (A.9) 

z — < oj r 
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